The Astrophysical Journal, 768:60 (1 lpp), 2013 May 1 
© 2013. The American Astronomical Society. All rights reserved. Printed in the U.S.A. 


doi:10.1088/0004-637X/768/l/60 


PHOTOIONIZATION MODELING OF OXYGEN K ABSORPTION IN THE INTERSTELLAR MEDIUM: 
THE CHANDRA GRATING SPECTRA OF XTE J1817-330 

E. Gatuzz 1 , J. Garcia 2,3 , C. Mendoza 1 ’ 4 , T. R. Kallman 3 , M. Witthoeft 3 , 

A. Lohfink 2 , M. A. Bautista 5 , P. Palmeri 6 , and P. Quinet 6 - 7 
1 Centro de Ffsica, Institute) Venezolano de Investigaciones Cientfficas (IVIC), P.O. Box 20632, 

Caracas 1020A, Venezuela; egatuzz@ivic.gob.ve, claudio@ivic.gob.ve 
- Department of Astronomy and Maryland Astronomy Center for Theory and Computation, University of Maryland, 

College Park, MD 20742, USA; javier@astro.umd.edu, alohfink@astro.umd.edu 
3 NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA; timothy.r.kallman@nasa.gov, michael.c.witthoeft@nasa.gov 
4 Centro Nacional de Calculo Cientlfico Universidad de Los Andes (CeCalCULA), Corporation Parque Tecnologico de Merida, Merida 5101, Venezuela 
5 Department of Physics, Western Michigan University, Kalamazoo, MI 49008, USA; manuel.bautista@wmich.edu 
6 Astrophysique et Spectroscopie, Universite de Mons-UMONS, B-7000 Mons, Belgium; palmeri@umons.ac.be, quinet@umons.ac.be 
7 IPNAS, Sart Tilman B15, Universite de Liege, B-4000 Liege, Belgium 
Received 2012 October 12; accepted 2013 March 6; published 2013 April 15 


ABSTRACT 

We present detailed analyses of oxygen K absorption in the interstellar medium (ISM) using four high-resolution 
Chandra spectra toward the X-ray low-mass binary XTE J 1 817-330. The 11-25 A broadband is described with 
a simple absorption model that takes into account the pile-up effect and results in an estimate of the hydrogen 
column density. The oxygen K-edge region (21-25 A) is fitted with the physical warmabs model, which is based 
on a photoionization model grid generated with thexSTAR code with the most up-to-date atomic database. This 
approach allows a benchmark of the atomic data which involves wavelength shifts of both the K lines and 
photoionization cross sections in order to fit the observed spectra accurately. As a result we obtain a column density 
of Ah = 1.38 ± 0.01 x 10 21 cirT 2 ; an ionization parameter of log £ = —2.70 ± 0.023; an oxygen abundance of 
A 0 — 0.6897q oVol an d ionization fractions of Oi/O = 0.91 1, On/O = 0.077, and Om/O = 0.012 that are in good 
agreement with results from previous studies. Since the oxygen abundance in warmabs is given relative to the solar 
standard of Grevesse & Sauval, a rescaling with the revision by Asplund et al. yields Ao = 0.952 4 q q 2 ®, a value 
close to solar that reinforces the new standard. We identify several atomic absorption lines — Key, K/3, and K y in O i 
and O n and Ka in O in, O vi, and O vn — the last two probably residing in the neighborhood of the source rather 
than in the ISM. This is the first firm detection of oxygen K resonances with principal quantum numbers n > 2 
associated with ISM cold absorption. 

Key words: atomic data - atomic processes - ISM: general - stars: individual (XTE J1 8 17-330) - X-rays: 
binaries - X-rays: general 
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1. INTRODUCTION 

X-ray studies of interstellar absorption have the potential to 
provide information that is not accessible by other techniques. 
This includes the relative abundances of a wide range of ion 

stages of interstellar cosmic elements (e.g., oxygen) and the 

abundances relative to H and He of elements with atomic 
number Z ^ 6. X-ray absorption can also provide signatures 
of the binding of these elements in molecules or solids; the 
inner-shell electronic transitions are key diagnostics since the 
ionization state or chemical binding shifts the line energies 
by a predictable amount. The ongoing effort to explore and 
exploit this effect in interstellar studies involves accumulation of 

synthetic spectra for inner-shell absorption of ions and neutrals, 
based on both atomic calculation and laboratory experiment, in 
addition to testing them against available observed interstellar 
X-ray absorption spectra. 

As part of the effort to accumulate the required atomic data, 
the atomic database of the xstar photoionization modeling code 
(Bautista & Kallman 2001) has been systematically improved 
in the past 10 years in order to study the K lines and edges in 
the high-quality X-ray astronomical spectra obtained from the 

Chandra and XMM-Newton satellite-born observatories. This 
database includes K-vacancy levels, wavelengths, A values, 
radiative and Auger widths, and high-energy photoionization 


cross sections for complete isonuclear sequences with Z ^ 30 
(see Palmeri et al. 2012 and references therein). 

With this new database it has been possible, for instance, 
to estimate the model efficiency of iron K-line emission and 
absorption in terms of the ionization parameter and density 
(Kallman et al. 2004). It has been shown that the centroid of 
the Ka unresolved transition array, the K/3 energy, and the ratio 
of the K«i to Ka 2 components are useful diagnostics of the 
ionization parameter. It was also found that the many strongly 
damped resonances below the K-ionization thresholds lead to 
edge smearing, which can certainly hamper the astrophysical 
interpretation of the absorption features. A synthetic spectrum 
has been modeled by Kallman et al. (2009) for the black hole 
X-ray transient GRO J 165 5 -40 which constrains the odd Z 
element abundances (11 ^ Z ^ 27) and outflow parameters 
of the associated warm absorber. By analyzing the oxygen 
K-absorption structure in the interstellar medium (ISM) from the 
XMM-Newton spectrum of the low-mass X-ray binary Sco X- 1 , 
Garcia et al. (2011) reproduced both the K edge and Oi Ka 
absorption line, thus evaluating the impact of the atomic data 
on the interpretation of observations. In the present report 
we continue this work with a more stringent benchmark of 
the xstar atomic database, made possible by Chandra spectra 
of a low-mass binary displaying imprints of the ISM oxygen K 
absorption with higher resolution. 
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Juett et al. (2004) measured high-resolution spectra of the 
interstellar oxygen K shell absorption edge from seven X-ray 
binaries— Cygnus X-l, Cygnus X-2, 4U 1636-53, 4U 1735-44, 
GX 9+9, 4U 1543-624, and 4U 1820-30— using the Chandra 
High Energy Transmission Grating Spectrometer (HETGS). 
The oxygen column density No was therein calculated using 
the optical depth at 21.7 A and the cross section at this 
wavelength assuming the oxygen abundance from Wilms et al. 
(2000). With this value, the hydrogen column density Ah 
was then determined. The \s-2p transitions of On and Om 
were identified, finding ISM abundances for Oil and Om 
relative to Oi of RiO.l and ^0.1, respectively. This became the 
first estimate of the Oii/Oi and Om/Oi abundance ratios in 
the ISM. 

By observing the oxygen K-shell edge toward galaxy clus- 
ters, Baumgartner & Mushotzky (2006) showed that the ISM 
hydrogen X-ray column density was in close agreement with 
the 21 cm radio value for columns less than approximately 
0.5 x 10 21 cm~ 2 . For higher column densities, estimates by 
Baumgartner & Mushotzky were higher by as much as a factor 
of 2.5, indicating substantial absorption besides that due to neu- 
tral hydrogen (probably from clouds of molecular hydrogen). 
An average ISM oxygen abundance of 0.99 solar was found rel- 
ative to the photospheric value (8.66 ± 0.05) by Asplund et al. 
(2004) which implied a high gas-to-dust ratio. 

Pinto et al. (2010) reported a spectral model suggesting that, 
along the line of sight toward an X-ray source (GS 1826- 
238), the ISM was more complex than just simple neutral 
gas. Their gas model consisted of three components: cold gas 
( T ^ 5.8-10 x 10 3 K), warm ionized gas ( T 1-6 x 10 4 K), 

and hot ionized gas (f w2x 10 6 K). It was also shown that the 
column densities Ah of these three components span two orders 
of magnitude, where the cold gas contributed ~90%-95%, the 
warm gas ~5%-10%, and the hot gas ~1% to the total A^ 1 . 
Costantini et al. (2012) analyzed the ISM using the bright X-ray 
binary 4U 1820-30, concentrating on the mildly ionized cold 
gas where the model consisted of mildly ionized gas, dust, and 
mildly ionized outflowing gas. In the oxygen edge region, dust 
compounds were included which did not show features below 
23.7 A. Finally, oxygen was found to be overabundant by a factor 
of 1 .23 solar. 

In order to model the oxygen K-absorption features in the ISM 
withxSTAR, we need a Chandra spectrum with high counts. In 
the catalog of low-mass X-ray binaries by Liu et al. (2007), 
we found in XTE J 1 8 17-330 a suitable source with an average 
count number of ~8 x 10 6 counts s -1 . Discovered by the Rossi 
X-Ray Timing Experiment (Remillard et al. 2006), XTE J1817- 
330 is a low-mass X-ray binary displaying a primary component 
which has an X-ray spectrum and variability behavior consistent 
with its identification as a black hole (Sala et al. 2007). Its 
location, / = 359.81 and b = —7.99 in Galactic coordinates, 
implies a line of sight in the Galactic center direction, a 
region that dominates the total mass of the Galaxy (Stasinska 
et al. 2012) providing a good opportunity to study the ISM 
absorption features. This source has been observed in different 
spectral ranges, namely the optical (Torres et al. 2006), near- 
infrared (D’Avanzo et al. 2006), radio (Rupen et al. 2006), 
and X-ray (Roy et al. 2011). Sala et al. (2007) fitted its X- 
ray spectrum finding an absorbing hydrogen column of 1.55 x 
10 21 cnL 2 and identifying the following oxygen absorption 
lines: Oi Ka(23.52A), Oil Ka(23.35A), Om Ka(23.13A), 
Oi K/3(22.91 A), and Ovu K/K21 .60 A). Gierlinski et al. 
(2009) estimated a hydrogen column density of A H = 1 .08 ± 


Gatuzz et al. 

0.05 x 10 21 cm -2 , and Roy et al. (2011) investigated the 
timing characteristics in the spectra using a column of Ah = 
1.2 xAo 21 cm~ 2 . 

The outline of the present paper is as follows. In Section 2, 
we describe thexSTAR warmabs physical model which is used 
to obtain the final spectrum fits, in particular the details con- 
cerning its atomic database. The data reduction aspects of the 
observations, e.g., the pile-up effect, are reviewed in Section 3, 
while the results of the broadband fits and resolved absorption 
features are given in Sections 4 and 5. Section 6 is dedicated 
to comparisons of the hydrogen column density that allow an 
estimate of its reliability. Section 7 describes the oxygen col- 
umn densities derived from our fits. Finally, the discussion and 
conclusions are presented in Section 8. 

2. XSTAR / WARM AB S 

xstar 8 9 is a code package designed to determine the physical 
conditions in partially ionized gases. It calculates the temper- 
ature, ionization, excitation, absorption, and emission due to 
neutral and all ionized species of elements with atomic num- 
ber Z 30. As many relevant physical processes as possi- 
ble are taken into account, assuming a time stationary balance 
among them and a Maxwellian electron velocity distribution. 
Important processes for many problems relevant to X-ray as- 
tronomy include photoionization, electron impact collisional 
ionization and excitation, and radiative and dielectronic recom- 
bination. More complete descriptions of xstar and its atomic 
database are given by Kallman & Bautista (2001) and Bautista & 
Kallman (2001). In this work we employ the warmabs utility 
which allows for the calculation of xstar models within the 
widely used X-ray spectral fitting package xspec. 9 It utilizes 
precalculated tabulations of the ion fractions and level popula- 
tions in order to calculate synthetic emission and transmission 
spectra. It assumes that the gas responsible for emission or ab- 
sorption has a uniform ionization and temperature throughout, 
although it is possible to superimpose two or more components 
and simulate a non-uniform situation. 

2.1. Atomic Data 

The K photoabsorption data sets for oxygen ions with 
electron numbers 2 ^ A ^ 8 contained in thexSTAR atomic 
database were generated by Garcia et al. (2005). Energies for 
both valence and K-vacancy fine-structure levels, transition 
wavelengths, gf values, and radiative and Auger widths were 
computed with the atomic structure packages autostructure 
(Eissner et al. 1974; Badnell 1986, 1997) andHFR (Cowan 
1981). Electron correlation effects were taken into account 
by means of configuration-interaction expansions within the 
n ^ 3 complexes, and relativistic corrections were included 
with a Breit-Pauli Hamiltonian. The accuracy of the theoretical 
level energies and wavelengths does not generally match that 
attained in experiments (see Table 1), and the scanty availability 
of measurements for the K shell of O ions with A > 3 limits 
the implementation of term-energy corrections and wavelength 
shifts that can be introduced in order to compensate for the 
theoretical shortcomings. Furthermore, as also shown in Table 1 , 
experimental misassignments, such as that quoted by Gu et al. 
(2005) for O vi, are not uncommon. (The discrepancy in O hi is 
believed to be a typo.) 


8 http://heasarc.nasa.gov/lheasoft/xstar/xstar.html 

9 http://heasarc.gsfc.nasa.gov/xanadu/xspec/ 
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Table 1 

Comparison of EBIT and Theoretical Wavelengths (A) for K Lines 


Ion 

Lower Level (J) 

Upper Level ( J ') 

EBIT 

Theory 

Ovi 

\s 2 2s(\/2) 

ls2s2p(l/2, 3/2) 

22.0194 ±0.0016 a 
22.374 ± 0.008 b 

22.00 c 

22.05 d 

22.03 e 

Ov 

\s 2 2s 2 (0) 

U2s 2 2p(l) 

22.374 ± 0.003 a 
22.370 ±0.010 b 

22.33 c 

22.35 d 

22.37 e 

Oiv 

\s 2 2s 2 2p(\/2, 3/2) 

ls2s 2 2p 2 (l/2, 3/2) 

22.741 ± 0.005 b 

22.78 c 

22.73 d 

22.15 s 

O hi 

\s 2 2s 2 2p 2 (\, 2) 

U2i 2 2p 3 (l) 

22.071 ±0.006 b 

23.08 c 

23.05 d 

23.07 e 


Notes. 

a EBIT measurement (Schmidt et al. 2004). 
b EBIT measurement (Gu et al. 2005). 
c hullac calculation (Behar & Netzer 2002). 
d /{-matrix calculation (Pradhan et al. 2003). 
e hfr calculation (Garcia et al. 2005). 


Table 2 

Chandra Observations Used in This Paper 


ObsID 

Date 

Exposure 

(ks) 

Instrument 

Grating 

Label 

6615 

2006 Feb 13 

18 

HETG-ACIS 

MEG 

Observation 1 

6616 

2006 Feb 24 

29 

HETG-ACIS 

MEG 

Observation 2 

6617 

2006 Mar 15 

47 

HETG-ACIS 

MEG 

Observation 3 

6618 

2006 May 22 

51 

HETG-ACIS 

MEG 

Observation 4 


High-energy photoionization cross sections in the K-edge 
region for oxygen ions were obtained by Garcia et al. (2005) 
withBPRM, the Breit-Pauli /f- matrix suite of codes (Scott & 
Burke 1980; Scott & Taylor 1982;Berringtonetal. 1987; Seaton 
1987). A conspicuous effect in the resonance structure of these 
cross sections that must be treated in detail is Auger damping. 
It is the result of the dominant spectator KLL Auger decay 
channels of a photoexcited K-vacancy state, say, 


— » 2p‘ M ~ 2 np + e~ 

0) 

— > [2s]2p ll ~ l np + e~ 

(2) 

-4- [2s] 1 2p ,1 np + e~ 

(3) 


over the participator Auger KLn channels 

[\s]2p^np — > 2/? /i_1 + e~ (4) 

-* [2s]2p lx + e~ , (5) 

which leads to edge smearing by resonances with symmetric 
profiles of nearly constant width. They are handled within the 
/f- matrix package by means of an optical potential (Gorczyca 
& Robicheaux 1999; Gorczyca & McLaughlin 2000). Again, 
the accuracy of the resonance energy positions is limited by 
the absence of measurements that would enable empirical 
adjustments of the many K-vacancy target thresholds (i.e., the 
series limits) in the close-coupling expansion. 



Figure 1. Light curves of the four observations, ObsIDs 6615, 6616, 6617, and 
6618, in units of counts s -1 . The approximately constant average counts s -1 
that prevails is an indication of a low degree of variability. 

(A color version of this figure is available in the online journal.) 

3. OBSERVATIONS AND DATA REDUCTION 

In Table 2 we list the specifications of the four observations 
of XTE J 1 8 17-330 obtained by Chcindra using the HETGS in 
combination with the Advanced CCD Imaging Spectrometer 
(ACIS). The oxygen edge is accessible with the HETGS 
exclusively via the Medium Energy Gratings (MEG). Given 
the large flux of this source (F = 1350 /rJ, integrated in the 
2-10 keV energy range), these observations were taken in 
continuous clocking mode which significantly increases the 
temporal resolution in order to minimize the pile-up effect 
(Cackett et al. 2008). Pile-up is an inherent feature of CCD 
detectors, such as those of the ACIS instrument that occurs 
when two or more photons are detected as a single event possibly 
causing a deformation in the level and shape of the continuum 
(Miller et al. 2006b). It is stronger in the MEG than in the 
HEG due to the lower dispersion and higher effective area of 
the former. 10 Although pile-up is usually not an issue near the 
oxygen edge, it may be for the continuum which is used to 
establish the hydrogen column density by fitting the 1 1-25 A 
interval; hence, we consider the pile-up effect in this range 
by applying the correction model simple_gpile2 . si (Hanke 
et al. 2009) before spectral fitting. Among the most popular 
computer packages for fitting X-ray data are xspec and isis, and 
we have used the latter (version 1.6.2 11 ) to include the pile-up 
model. 

Even though the spatial resolution is reduced to one dimen- 
sion when the continuous clocking mode is used, there is essen- 
tially no difference with the timed exposure mode in the back- 
ground extraction or in the data reduction process. Hence, we 
have reduced the data sets using the standard CIAO threads. 12 
In some cases, the zero-order data were not telemetered; for 
these observations the zero-order position was then estimated 
by finding the intersection of the grating arms using the f indzo 
algorithm. 13 

Figure 1 shows the light curves of the four observations 
ObsID 6615, 6616, 6617, and 6618 in units of counts s _1 . 


10 http://cxc.harvard.edu/ciao/ahelp/acis_pileup.html 

1 1 http://space.mit.edu/cxc/isis/ 

12 http://cxc.harvard.edU/ciao4.4/threads/gspec.html 

13 http://space.mit.edu/cxc/analysis/findzo/ 
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Figure 2. Chandra MEG spectra of the X-ray binary XTE J1 8 17-330 simultaneously fitted in the 1 1-25 A region using the TBnew model, (a) ObsID 6615. (b) ObsID 
6616. (c) ObsID 6617. (dl ObsID 6618. See Table 2 for further details pertaining to these observations. In each panel, the black data points are the observations, while 
the solid red lines correspond to the best-fit models. Residuals are shown as histograms for each case in units of x 2 - Note the large residuals in the oxygen K-edge 
region, particularly in the neighborhood of the O I and O n Key lines at ^23.50 A and ^23.35 A, respectively. 


(A color version of this figure is available in the online journal. 1 


It can be seen that the average counts s _l is approximately 
constant for all, and therefore we use the averaged spectra for 
this work. We have carried out the following spectral fitting 
procedure: first, each spectrum was rebinned with a 0.1 A bin 
in order to fit the continuum in the 1 1-25 A wavelength range 
so as to find the hydrogen column density using the TBnew 
model; then, the default spectral resolution was restored, and the 
oxygen edge region (21-24 A) was analyzed with the warmabs 
model freezing the column density at the value obtained in the 
preceding step. 

4. BROADBAND FIT 

Following Yao & Wang (2006), we rebinned each spectrum 
with a 0. 1 A bin size throughout the spectral range to decrease 
the number of absorption lines and improve the broadband 
fit. We fitted the four observations simultaneously in the in- 
terval 1 1-25 A using the simplegpile2 (TBnew (powerlaw)) 
model, where simplegpile2 is a convolution model to account 
for pile-up effects and TBnew is an X-ray absorption model that 
includes abundances for elements from H to Ni (Wilms et al. 
2000). Also we added Gaussians to fit the remaining absorption 
lines after rebinning. We used abundances specified in Wilms 
et al. (2000) and cross sections by Verner et al. (1996). The 


O, Ne, and Fe abundances were handled as free parameters, 
and for the analysis x 2 statistics were employed. We fixed the 
absorption parameters in all the observations, and varied the 
power-law, pile-up, and Gaussian parameters for each case. 

Figure 2 shows the fit in the interval 1 1-25 A using the 
simplegpile2 (TBnew (powerlaw)) model. Although the fit 
is carried out simultaneously for the four observations, they are 
plotted separately to improve clarity. Panels (a), (b), (c), and 
(d) correspond to observations ObsID 6615, 6616, 6617, and 
6618, respectively (see Table 2). In each panel, black data points 
correspond to the observation, while the solid red lines represent 
the best-fit model for each case, and the base plots show the 
fit residuals in units of / 2 . In each spectrum the K edges of 
Ne and O are, respectively, located at «sl4.3A and ^23.2 A. 
The Oi and Oil ls-2p absorption lines at 23.5 A and 23.35 A 
are also clearly observed. We improved the fit by including 
Gaussians to model the outstanding absorption lines: Neix Ka 
(«13.43 A), Ne m K a («G4.50 A), Ne ii K a («14.60 A), O viii 
Lya (^18.95 A), O vn Ka (^21.58 A), and O ii Ka (^23.34 A). 
Residuals indicate that the model fits the data satisfactorily with 
deviations near the Ne edge («sl4.3 A) due to background (i.e., 
clocking the whole chip) or an unaccounted for hot pixel. We 
also find high residuals in the oxygen absorption region; in 
particular, the O i K« and O ii K« are not well modeled. 
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Table 3 

Broadband Simultaneous Fit Parameters 


Model 

Parameter 

Value 

TBNew 

Ah (10 21 cm- 2 ) 

1 66 +003 
1 • oo -0.04 

TBNew 

Neon abundance (An c ) 

2.258 

O 

cc 

o 

o 

-H 

TBNew 

Iron abundance (Ap e ) 

1.302 

± 0.050 

TBNew 

Oxygen abundance (Ao) 

1.178 

it 0.022 

Power law 

Normalization 

9.416 

± 0.064 a 



7.039 

±0.041 b 



4.613 

± 0.042 c 



1.481 

±0.001 d 

Power law 

Photon index 

1.545 

± 0.030 a 



1.618 

± 0.030 b 



1.735 

± 0.029 c 



2.242 

oo 

co 

o 

o 

-H 

SimpleGpile2 

f> 

0.050 

±0.001 a 



0.051 

±0.001 b 



0.050 

±0.001 c 



0.049 

±0.001 d 

Reduced chi-square 

x 2 

1 

.208 


Notes. Abundances relative to the solar values of Wilms et al. (2000). 
a ObsID 6615. 
b ObsID 6616. 
c ObsID 6617. 
d ObsID 6618. 

Fit results for the individual free abundances, column density, 
power law, and pile-up are presented in Table 3. The fit statistics 
give a reduced chi-square of / 2 = 1.208. Differences in the 
power-law parameters are due to the count numbers in each 
observation. We obtain abundances for O, Ne, and Fe higher 
than solar (Wilms et al. 2000), especially for Ne (a factor of 
two), and a column density of Nu — 1 .66^ q 3 4 x 10 21 cm -2 
which is somewhat larger than those estimated in other spectral 
fits and by measurements of the 21 cm line. This will be further 
discussed in Section 6. 

The nonlinear pile-up convolution model simple_gpile2 . si 
exponentially reduces the predicted count rate Rj(X) (in units 
of counts s^ 1 A” 1 ) according to 

R(k) = R(k)exp(-pR w (k)), (6) 

where R to t (X) is the total spectral count rate and /3 is treated as a 
fit parameter. The model is only applicable to first-order grating 
spectra, although it employs information from the second- and 
third-order grating spectra. 14 With this model we obtained an 
average j3 = 0.050 according to which the spectra have a pile- 
up degree greater than 25% for 14 ^ A. ^ 16 A with the highest 
value of p — 42% at ~ 1 4.3 A for the observation with the 
largest number of counts (ObsID 6615). This value shows the 
need to include the pile-up effect in the 1 1-25 A region to obtain 
a good fit. However, for all the observations, we have a pile-up 
degree lower than 5% in the 21-24 A wavelength range which 
corresponds to the oxygen absorption region and, therefore, its 
effects were ignored. 

5. ABSORPTION FEATURES 

In order to identify all the absorption features present in the 
21-25 A wavelength region, we first applied a functional model 
consisting of a power law for the continuum and Gaussian pro- 
files to describe the absorption lines. Figure 3 shows the spectral 


14 http://cxc.harvard.edu/ciao/ahelp/acis_pileup.html 
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Table 4 

Absorption Line Assignments 


Ion 

Transition 

A a 

(A) 

k b 

(A) 

O vn 

Is 2 1 5-[l.v]2p l P° 

21.589 ±0.003 

21.60 ±0.06 

Ovi 

2s 2 S-[ls]2s2p 2 P° 

22.022 ± 0.003 


On 

2 p 3 4 S°~[ls]2p 3 4p 4 P 

22.101 ±0.005 


On 

2p 3 4 S°-[ls]2p 3 3p 4 P 

22.280 ± 0.003 


Oi 

2 p 4 3 P-[ls]2p 4 ( 2 P)4p 3 P° 

22.609 ± 0.007 


Oi 

2 p 4 3 P-[ls]2p 4 ( 2 P)3p 3 P° 

22.686 ± 0.004 


Oi 

2 p 4 3 P-[ls]2p 4 ( 4 P)4p 3 P° 

22.790 ± 0.001 


Oi 

2p 4 3 P-[ls]2p 4 ( 4 P)3p 3 P° 

22.884 ± 0.004 

22.91 ±0.03 

Oiii 

2 p 2 3 P-[li]2p 3 3 P° 

23.054 ± 0.001 

23.13 ±0.09 

Oiii 

2 p 2 3 P-[li]2p 3 3 D°, 3 S° 

23.104 ±0.005 

23.13 ±0.09 

On 

2 p 3 4 S°-[ls]2p 44 P 

23.358 ± 0.002 

23.35 ± 0.03 

Oi 

2 p 4 3 P-[li]2p 5 3 P° 

23.502 ±0.001 

23.52 ± 0.02 


Notes. 

a Present work. 
b Sala et al. (2007). 


fit of the four Chandra MEG observations of XTE J 18 17-330 in 
the oxygen absorption region. Although shown in separate pan- 
els for clarity, the fit was performed simultaneously for the four 
observations. Nevertheless, the power-law parameters (photon 
index and normalization) are allowed to vary independently for 
each data set. Several absorption features are clearly observed, 
the most outstanding being the O i and O ii Key lines. The base 
plot in each panel shows, in units of / 2 , the fit residuals as a his- 
togram. High residuals are found around 24.8 A, and the model 
underestimates the number of photons around 23.2 A. 

Table 4 shows a list of the lines obtained with the func- 
tional fit including the centroid wavelengths derived from the 
Gaussian profile parameters. It also lists the wavelengths for 
the lines measured by Sala et al. (2007). Note that our wave- 
lengths are within the error bars of the values reported by these 
authors, which are much larger given the lower spectral reso- 
lution of the XMM-Newton instruments. We obtain a reduced 
chi-square of / 2 = 1.434. The Oi Key is found at 23.502 ± 
0.001 A which has been previously declared to be of ISM origin 
(Sala et al. 2007). Furthermore, we have also detected four other 
Oi absorption lines corresponding to resonances with higher 
principal quantum number n: two pairs of Kf J > (n = 3) and Ky 
(n = 4) resonances associated with the 4 P and 2 P K-hole core 
states. These are found at 22.884 ± 0.004 A, 22.790 ± 0.001 A, 
22.686 ± 0.004 A, and 22.609 ± 0.007 A, respectively. In this 
set, only the first Kp (belonging to the 4 P core state) has been 
previously detected at 22.91 ± 0.03 A in the XMM-Newton ob- 
servation analyzed by Sala et al. (2007). This is then the first 
firm detection of oxygen K resonances with n > 2 associated 
with ISM cold absorption. 

The line at 23.358 ± 0.002 A is also of ISM origin (Juett 
et al. 2004) and corresponds to On Ka. As in the case of Oi, 
resonances with n > 2 were also detected for O n. These were 
found at 22.280 ±0.003 A and 22.101 ±0.005 A corresponding 
to Kfi and Ky , respectively. Although the O in Key resonance 
is a triplet, only two absorption features were detected at 
23.104 ± 0.005 A and 23.054 ± 0.001 A __ while Sala et al. 
(2007) only found one line at 23.13 ± 0.09 A identified as Om 
Key . Finally, we have also located Key transitions for both O vi 
(22.022 ± 0.003 A) and O vn (21.589 ± 0.003 A), from which 
only the strongest of the two (O vn) was reported by Sala et al. 
(2007) at 21.609 ±0.06 A. 
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Figure 3. Spectral fit of the Chandra MEG observations of XTE J 1 8 17-330 in the oxygen absorption region (21-25 A) using a simple functional model (power law 
and several Gaussian profiles), (a) ObsID 6615. (b) ObsID 6616. (c) ObsID 6617. (d) ObsID 6618. The observed absorption features are interpreted as the neutral 
O K-edge (~23.1 A); the Key, K ft, and K y lines of Oi («s23.50A. ^22.88 A, and ^22.68 A, respectively); the Key, K/3, and K y lines of On (^23.35 A, ^22.28 A. 
and 5=»22. 10 A, respectively); and the Key lines of O in, O vi, and O vii (~23. 1 0 A, RJ22.02 A, and ~21 .58 A, respectively). See Table 4 for more details. 

(A color version of this figure is available in the online journal.) 


The functional fit described above illustrates the complexity 
of oxygen K absorption. There are (at least) two different 
regimes: a cold component of mostly neutral gas and a hot 
component with highly ionized gas. Previous studies have 
suggested a multiple-phase ISM; however, it is not clear if 
the two high-ionization lines from O vi and O vii are intrinsic 
to the source or in fact arise in the ISM. Thus we direct our 
study to the cold phase. Furthermore, since cold absorption 
is not exclusively due to neutral oxygen, this means that a 
model such as TBnew will not completely represent the observed 
features. We have then simultaneously fitted the four exposures 
in the 21-25 A wavelength range using the powerlaw*warmabs 
physical model. The advantage of using warmabs is that it 
contains the most recent atomic data for all the ions in the 
oxygen isonuclear sequence, thus enabling the fit of all the 
absorption lines from O i, O ii, and O ill. While implementing 
this model, the column density is held fixed at the value obtained 
in the broadband simultaneous fit (Ah = 1.66 x 10 21 cm -2 ). 
The ionization parameter and the oxygen abundance are treated 
as free parameters but the same for all observations, while 
the abundances of all the other chemical elements are held 
fixed at their solar values. The quantities describing the power 
law (photon index and normalization) are also taken as free 
parameters but independently for each observation. 


Figure 4 shows the results of this fit. As before, we present 
for clarity the spectra with the best-fit model in separate panels, 
although the fit is performed simultaneously. Since most of the 
relevant absorption features in the spectra are due to neutral, 
singly, and doubly ionized O ions, the ionization parameter that 
describes the fit tends to be somewhat low: log f = —2.729. In 
this case, the warmabs model does not show any absorption due 
to O vii; therefore, this line is still represented with a Gaussian 
profile. We have found large residuals around most of the oxygen 
absorption features, the most prominent being around the O i and 
Oil Ka lines. The fit statistics give a poor reduced chi-square 
of x 2 — 1-769 which forced us to revise the atomic database. 

A close examination of the resonance positions from the 
photoabsorption cross sections used in warmabs in comparison 
with those obtained from the spectral fits with Gaussian profiles 
reveals the inaccuracies of the atomic data. Figure 5(a) shows 
with solid colored lines the atomic cross sections used in 
warmabs for Oi, Oil, and Oiii (Garcia et al. 2005). Vertical 
dashed lines are placed at the positions of the absorption features 
found with the functional fit described above (see Table 4). A 
comparison shows that not only the positions of the Oi and 
O ii Ka resonances are displaced with respect to the observed 
data but also those with n > 2. The solid black line represents 
the experimental Oi photoabsorption cross section measured 
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Figure 4. Spectral fit of the Chandra MEG observations of XTE J1 8 17-330 in the oxygen absoiption region (21-25 A) using a powerlaw*warmabs physical model, 
(a) ObsID 6615. (b) ObsID 6616. (c) ObsID 6617. (d) ObsID 6618. 

(A color version of this figure is available in the online journal.) 
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Figure 5. Comparison of the theoretical photoabsorption cross sections for O I (red), O n (green), and O hi (blue) computed by Garcia et al. (2005) that are implemented 
in the warmabs model. The black solid line is the laboratory measurement by Stolte et al. (1997) showing peaks of Oi and 02- The vertical dashed lines are placed at 
the wavelengths of the observed absorption lines (see Table 4). Panel (a) displays the original cross sections, while panel (b) shows the same curves after the wavelength 
shifts are applied to both O I and O n (see the text for details). 

(A color version of this figure is available in the online journal.) 


by Stolte et al. (1997; unfortunately, there are no laboratory 
measurements available for the other species) which displays a 
similar wavelength shift with respect to the observed resonance 
positions. This wavelength offset has been previously discussed 
by Juett et al. (2004) based on the K« position observed 
in the Chandra spectra from several sources, all of which 


match the position we have found using Gaussian profiles. For 
Oi, the comparison of the theoretical cross section with the 
laboratory measurement shows that, except for the K« peak, 
all the other resonances agree very well. This implies that the 
relative position of the Ka with respect to the other resonances 
also needs to be adjusted. This is also the case for Oil since 
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Figure 6. Spectral fit of the Chandra MEG observations of XTE J1 8 17-330 in the oxygen absorption region (21-25 A) using a powerlaw*warmabs physical model, 
after the wavelength corrections in the Oi and On atomic cross sections are introduced, (a) ObsID 6615. (b) ObsID 6616. (c) ObsID 6617. (d) ObsID 6618. 

(A color version of this figure is available in the online journal.) 


the wavelength difference of the theoretical and observed Ka 
positions is not the same as those for the rest of the resonances 
of the same ion. Consequently, we have decided on the one hand 
to adjust the positions of the Ka resonances (which are treated 
as lines in warmabs), and on the other hand, to shift all the 
cross sections for both species so as to obtain the best possible 
agreement with the observed lines. The new positions for the O i 
and Oil Ka lines are 23.502 A and 23.343 A, while the cross 
sections were shifted toward shorter wavelengths by 0.033 A 
and 0.079 A, respectively. No correction was applied to the O hi 
curve since the observed lines are weaker, in addition to the 
fact that in this case the Ka resonance is a triplet. Figure 5(b) 
shows the atomic absorption cross sections after the wavelength 
corrections are applied. In other words, the Chandra observed 
line positions suggest not only that the theoretical resonance 
positions from Garcia et al. (2005) must be adjusted, but also 
that the wavelength scale of the only available experiment for 
O i (Stolte et al. 1997) needs to be shifted by 33 mA. 

Possible effects due to the instrument wavelength calibration, 
model uncertainties, and Doppler shifts due to motion of the 
gas in the ISM are important for accurate line identification. 
The instrumental resolution for the first-order spectra from 
MEG is AA ~ 23 mA, which is better than the largest shift 
required for our theoretical cross sections (33 mA). In fact, 
according to Garcia et al. (2005), the uncertainties in the 
theoretical cross sections for O ions with electron occupancies 
N ^ 4 were estimated to be ~50 mA, comparable with the 


wavelength uncertainties in the laboratory measurements of 
Stolte et al. (1997). Assuming a velocity dispersion for the ISM 
of ^200 km s -1 (luett et al. 2004), the largest wavelength shift 
due to Doppler effects for the O i Ka is AA = 14 mA, which is 
smaller than the instrumental resolution. Hence, we can safely 
ignore the effects of ISM gas motion in the line of sight. 

Figure 6 shows the fit results with warmabs after the theoret- 
ical line positions have been corrected, where the residuals are 
now significantly reduced and more evenly distributed. This is 
reflected in better fit statistics (reduced / 2 = 1.245) and a 
well-constrained ionization parameter (log £ = —2.699 ± 
0.023) which is consistent with the low-ionization oxygen lines 
observed. As before, the Ovu Ka resonance at 21.589 A is 
not well represented with this model and must be fitted with a 
Gaussian profile. We have also noticed that the absorption fea- 
ture at 22.022 A is not entirely due to the Om KA resonance 
and is thus fitted with a Gaussian profile, confirming a blend 
of both Oiii KA and Ovi Ka. In this final analysis, we also 
treated the column density value as a free parameter in order 
to increase the accuracy of the fit. We find a column density 
of /V H = 1.38 ± 0.1 x 10 21 cm -2 that is lower than the value 
previously mentioned in connection with the broadband fit (see 
Section 4). The oxygen abundance obtained with the warmabs 
model is A 0 = 0.689+oq 1 1 5 0 relative to solar, also lower than 
the value for the best fit with TBnew (Ao = 1.178 ± 0.22). 
All the best-fit parameters of the warmabs are summarized in 
Table 5. 
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Table 5 

warmabs Simultaneous Fit Parameters 


Model 

Parameter 

Value 

Warmabs 

N h (10 21 cm- 2 ) 

1.38 ±0.01 

Warmabs 

Log ionization parameter (log£) 

-2.699 ± 0.023 

Warmabs 

Oxygen abundance (Ao) 

n Aso + t).0l5 

0.689_ 0010 

Power law 

Normalization 

2.384 a 

2.544 b 



1.51 l c 
0.685 d 

Power law 

Photon index 

3.044 a 

2.479 b 



2.720 c 

2.617 d 

Reduced chi-square 

x 2 

1.245 


Notes. Oxygen abundance relative to the solar value of Grevesse & Sauval 
(1998). 

a ObsID 6615. 
b ObsID 6616. 
c ObsID 6617. 
d ObsID 6618. 

Table 6 

Hydrogen Column Density Comparison 


Method A , n(xl0 21 cm 2 ) 


TBnew fit a 

1 66 +003 
1,0 — 0.04 

warmabs fit 3 

1.38 ± 0.1 

21 cm survey b 

1.58 

2 1 cm survey c 

1.39 

Spectral model fit d 

1.55 ± 0.05 

Spectral model fit e 

1.2 

Spectral model fit f 

0.88-0.97 

Spectral model fit g 

1. 0-3.0 

Spectral model fit h 

1.08 ± 0.05 


Notes. 

a Present work. 
b Dickey & Lockman (1990). 
c Kalberla et al. (2005). 
d Sala et al. (2007). 

e Rykoff et al. (2007), Gierliriski et al. (2008), Roy et al. (2011). 
f Miller et al. (2006a). 

8 Torres et al. (2006). 
b Gierliriski et al. (2009). 

6 . HYDROGEN COLUMN DENSITY 

Table 6 shows a comparison of the hydrogen column den- 
sity obtained in this work with previous estimates including 
measurements using the 21 cm line. By means of the simulta- 
neous broadband fit described in Section 4, we have obtained 
a column density of Ah = 1 . 66^004 x 10 21 cm 2 which is 
somewhat larger than those quoted in previous spectral fits. 
Miller et al. (2006a) reported a hydrogen column density of 
Ah = 8 . 8-9.7 x 10 20 cm ' 2 using a 50 ks Chandra spec- 
trum, a value smaller than those obtained from Galactic Hi 
surveys, e.g., 1.58 x 10 21 cm -2 (Dickey & Lockman 1990) and 
1.39 x 10 21 cnT 2 (Kalberla et al. 2005). Using the equivalent 
widths (EWs) of several interstellar bands in the optical spec- 
trum, Torres et al. (2006) found the column density to be in 
the range of /V H = 1-3 x 10 21 cm 2 . Moreover, both Rykoff 
et al. (2007) and Gierliriski et al. (2008) gathered acceptable fits 
of several Swift X-Ray Telescope (XRT) spectra with a fixed 
Galactic column density of Ah = 1.2 x 10 21 cm' 2 . Gierliriski 
et al. (2009) used the same XRT spectra combined with data 
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Table 7 

Equivalent Width Comparison for the Oxygen Ka Transitions in the ISM 


Source 

0 1 Ka? 

OuKa 

O ill Ka 

Ovi Ka 

O vii Ka 


(mA) 

(mA) 

(mA) 

(mA) 

(mA) 

XTE J1817-330 2 

51 ±5 

52 ±5 

16 ±5 

14 ±5 

54 ±4 

XTE J1817-330 b 

70 ± 20 

40 ± 20 

30 ± 30 


40 ± 30 

Cyg X-2 C 

67 ± 11 

26 ± 17 

9 ± 8 



Cyg X-2 d 

50 ±3 

28 ±6 


11 ± 2 


4U 1543-624 c 

68 ± 11 

41 ± 11 

18 ± 5 



4U 1820-30 c 

70 ±20 

41 ± 17 

20 ± 10 



4U 1735-44° 

80 ±30 

40 ± 20 

24 ± 10 



GX 9+9 c 

80 ±30 

50 ± 19 

27 ± 13 



4U 1636-53° 

90 ±30 

60 ± 20 

33 ± 13 



Cyg X-l (ObsID 3407)° 

95 ± 17 

20 ± 1 1 

35 ± 15 




Notes. 

3 Present work. 
b Sala et al. (2007). 
c Juett et al. (2004). 
d Yao et al. (2009). 

from UVOT (the Ultraviolet and Optical Telescope) to derive a 
column density of Ah = 1.08 ± 0.05 x 10 21 cm' 2 from broad- 
band spectral fits. 

The agreement between TBnew and Sala et al. (2007) is 
around 7% despite the fact that their fits include data from 
XMM-Newton EPIC-Pn (0.6-10.0 keV), RGS1 (0.3-2.0 keV), 
and OM covering a much wider energy range; furthermore, 
they model the foreground absorption with TBabs, an older 
version of TBnew (Wilms et al. 2000). On the other hand, 
the present fits of the oxygen K region (21-25 A) with the 
warmabs physical model yield a lower column density (Ah = 
1 .38 ± 0. 1 x 10 21 cm” 2 ) than that derived from our broadband 
fit and TBnew. This discrepancy is expected due to the different 
atomic data sets involved, particularly since TBnew only includes 
the photoabsorption cross section for neutral oxygen. Taking 
into account the scatter of the previous column densities in the 
literature and the fact that we have adequately modeled the 
absorption features from Oi, On, and Om, we would expect 
the value derived from warmabs to be more reliable than that 
from TBnew. 

7. OXYGEN COLUMN DENSITY 

In order to derive the O column density along the line of sight 
of XTE J 1 817-330, we have calculated the curve of growth for 
each of the oxygen Ka transitions using the atomic data from 
Garcia et al. (2005), and EWs have been obtained using the 
functional fit described in Section 5. Table 7 shows a comparison 
between the present Ka EWs for 0 1 , 0 11 , O in, O vi, and O vn 
and those by Sala et al. (2007) for XTE J1817-330. EWs listed 
in Juett et al. (2004) for seven other sources are also included 
as well as those by Yao et al. (2009) for Cygnus X-2. For 0 1 
Ka, our EW (51 ± 5 mA) is around the lower limits in Sala 
et al. (2007) and Juett et al. (2004), and for 0 11 Ka, the present 
value of 52 ± 5 mA is within the range that Juett et al. (2004) 
associated with the ISM. However, we find the Ka lines in both 
0 1 and 0 11 to be saturated which could influence our derived 
EWs. Furthermore, the Ka EWs in Om, Ovi, and Ovii are in 
good agreement with the values obtained by Juett et al. (2004), 
Sala et al. (2007), and Yao et al. (2009). 

Table 8 shows a comparison of the oxygen column densi- 
ties obtained from the EWs and the warmabs model fit. As 
previously mentioned in Section 5, for the ionization param- 
eter found in our best warmabs fit, only Oi, Oil, and Om 
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Table 8 

Oxygen Column Density Comparison 
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Source 

Oi 

(x 10 17 cm' 2 ) 

On 

(x 10 17 cm' 2 ) 

Oiii 

(xlO 17 cm' 2 ) 

Ovi 

(x 10 15 cm' 2 ) 

O vu 

(xlO 15 cm' 2 ) 

XTE J1817-330 a 

3.18 ±0.79 

2.18 ±0.54 

0.35 ±0.08 

5.06 ± 2.53 d 

45 ± 18 d 

XTE J1817-330 b 
PG 1302-102 c 
Mrk 1383 c 
ESO 141-G55 C 
PKS 2005-489 c 

5.85 ± 1.75 

0.49 ±0.12 

0.07 ± 0.02 

5.99 ± 1.9 e 
71.12 ±29.15 f 

0.16 ±0.03 
0.38 ±0.07 
0.31 ±0.06 
0.60 ± 0.01 

988 ± 405 e 
4484 ± 1838 f 


Notes. 

a Calculated from the EWs and the curves of growth (see Section 7). 

b Derivated from the warmabs model fit using the solar value of Grevesse & Sauval (1998, see Section 5). 
c Obtained from FUSE data by Savage et al. (2003). 
d Using a velocity dispersion v = 200 km s' 1 . 
e Using a velocity dispersion v = 100 km s' 1 . 
f Using a velocity dispersion v = 20 km s' 1 . 


are taken into account in the model. As a reference, we have 
also included the oxygen column density for Ovi from Sav- 
age et al. (2003), obtained from Far-Ultraviolet Spectroscopic 
Explorer (FUSE) data for four sources located near the same 
Galactic latitude as XTE J 1817-330. The curve-of-growth val- 
ues for Oi, On, and Om result in a total column density of 
No = (Voi + AW + AWi = 5.71 x 10 17 cm' 2 , in relatively 
good agreement with that derived from the ionic fractions in 
the warmabs fit, namely No = 6.41 x 10 17 cm' 2 . However, 
the ratios Oi/Oii are significantly different: 1.4 and 11.9, re- 
spectively. Due to line saturation, the measured EWs are likely 
to be underestimated, in particular for the Oi Ka line. This 
leads to a higher uncertainty in the column and ion fractions 
estimated from the EW measurements. On the other hand, since 
the parameters in the warmabs fit are constrained by all the 
lines as well as the K edge, we are more confident about their 
reliability. In the case of the O vi and O vn Ka transitions, be- 
cause the EWs are on the flat section of the curve of growth, they 
strongly depend on velocity dispersion. For O vi Ka, the oxygen 
column density spans from /Vo VI = 71.12 ± 29.15 x 10 15 cm' 2 
to AW = 5.06 ± 2.53 x 10 15 cm' 2 when using dispersion ve- 
locities of v — 20 km s' 1 and v = 200 km s' 1 , respectively. The 
smallest value is still about an order of magnitude larger than that 
reported by Savage et al. (2003), making it difficult to constrain 
the velocity dispersion. Also, no real improvement is achieved 
by considering larger dispersions, as the column density is 
hardly reduced once the velocities are larger than 200 km s' 1 . 
This may be an indication that the observed O vi line is likely to 
originate from the neighborhood of XTE J 1 8 17-330 rather than 
from the ISM. In the case of O vu, the oxygen column density 
varies from AWn = 4-48 ± 1-84 x 10 18 cm' 2 (v = 20 km s' 1 ) 
to AW, = 4.5 ± 1.8 x 10 16 cm' 2 ( v = 200 km s' 1 ). 

8. DISCUSSION AND CONCLUSIONS 

There are several interesting conclusions emerging from this 
work that configure a set of useful guidelines in the analysis 
of ISM absorption features in Chandra and XMM-Newton 
spectra. In order to resolve such spectral signatures, bright 
sources with high count rates are desirable. In the present 
XTE J1 817-330 study case, we were fortunate to find four 
good quality spectra that were fitted simultaneously so as to 
obtain improved statistics. However, it was realized that the 
pile-up effect must be taken into account in order to obtain 


reasonably reliable hydrogen column densities. This effect was 
treated adequately with the nonlinear pile-up convolution model 
simple_gpile2 . si, whereby the pile-up was estimated by 
fitting the continuum to be on average around 25%. The resulting 
hydrogen column density is Ah = 1. 66*004 x 10 21 cm' 2 , 
which is somewhat (<7%) higher than those from Kalberla 
et al. (2005) and Sala et al. (2007), the former obtained from 
21 cm data. 

The absorption features in the oxygen region (21-25 A) were 
first picked up with a functional model (power law and Gaussian 
profiles) which yielded Ka lines of Oi, Oil, Om, Ovi, and 
O vu previously observed by Sala et al. (2007), but for the first 
time, also K resonances of O i and O ii with principal quantum 
number n > 2. These line identifications at least confirm 
a two-phase plasma with both neutral (mostly) and highly 
ionized components. In the present work we have been mainly 
concerned with the former. In this respect, since neutral oxygen 
is not the sole charge state responsible for the cold oxygen 
absorption, the four exposures are then simultaneously fitted 
with a representative physical model referred to as warmabs, 
which contains the atomic data for the oxygen isonuclear 
sequence computed by Garcia et al. (2005). In this step, 
the hydrogen column density is held fixed at the previously 
determined value of Ah = 1.66 x 10 21 cm' 2 , while the 
ionization parameters and oxygen abundance are treated as 
free parameters. A detailed comparison of this model with 
observations brings forth the inaccuracies of the atomic data, 
which led us to introduce wavelength shifts in both the Ka line 
positions and photoabsorption cross sections of O i and O ii. 

The scope of this corrective procedure is severely limited by 
the general poor availability of laboratory data for the oxygen 
inner shell, but it certainly leads to improved fits and statistics, 
in particular, to a well-constrained ionization parameter of 
logf = —2.699 ± 0.023 consistent with the absorption lines 
observed. In order to further improve the fit, the column 
density is allowed to vary resulting in a more reliable value 
(Ah = 1.38±0.1 x 10 21 cm' 2 ) which is within 13% of previous 
estimates (Kalberla et al. 2005; Sala et al. 2007). Furthermore, 
the fit yields an oxygen abundance of Aq = 0.689 A ( q ' ll 1 m and 
ionization fractions of Oi/O = 0.911, Oii/O = 0.077, and 
Om/O = 0.012. Regarding the oxygen abundance, it may be 
clarified that the quoted value is given relative to the solar 
standard (8.83 ± 0.06) of Grevesse & Sauval (1998). If it is 
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rescaled to the recently revised standard by Asplund et al. (2009) 
of 8.69 ± 0.05, then an ISM abundance (A 0 = 0.952+° 0 0 0 2 1 0 3 ) 
close to solar is obtained which is in excellent accord with the 
averaged value by Baumgartner & Mushotzky (2006), and can 
be included as further support for the solar abundance revision. 
Moreover, our ionization fractions are in good agreement with 
those found by Juett et al. (2004) of O ii/O i & 0.1 and O m/O i 
< 0 . 1 . 

The functional fit allows us to derive oxygen column densities 
by means of a curve-of-growth method using the EWs for 
each of the observed Ka absorption lines, which can then be 
compared with the values obtained from the warmabs fit. Even 
though we find by both methods a relatively consistent total 
column density (5.71 x 10 17 cm -2 in the curve of growth and 
N — 6.41 x 10 17 cm~ 2 in the warmabs model fit), the Oi/Oii 
ratio for the Key lines is significantly lower (an order of 
magnitude) in the former. This is attributed to the high saturation 
level in both lines which leads to underestimates of the measured 
EWs, and we are therefore confident that the results from the 
warmabs model fit are more reliable. In the case of Ovi, the 
width of the line strongly depends on the velocity dispersion of 
the gas. For the velocity dispersion values expected in the ISM 
( v ~ 20-200 km s -1 ), the EWs obtained with our functional fit 
yield oxygen column densities far above from the FUSE values 
by Savage et al. (2003). This may be an indication that the O vi 
observed absorption in the XTE J18 17-330 spectra is bound to 
occur near the source rather than in the ISM. 

We find acceptable fits to the observed oxygen K-edge 
spectrum with a model consisting of atomic oxygen ions if 
the energies are slightly adjusted. However, it is also of great 
interest to detect or set limits on features due to other phases 
of interstellar oxygen and other elements, namely molecules 
or solids (dust). It is expected that these phases comprise a 
significant fraction of the total oxygen column in the ISM. 
Evidence for this includes observed depletions of atomic oxygen 
along many lines of sight relative to certain elements (Jenkins 
2009; Yao et al. 2009; Whittet 2010), together with observations 
of molecular features in other wavelength bands. An example 
of a potentially detectable signature of such material is the K 
resonance of the CO molecule expected at a rest wavelength 
of 23.21 A (Barrus et al. 1979). This line is separated from the 
nearest feature we clearly detect by approximately 100 mA, 
which can be compared with the magnitude of the wavelength 
shifts we adopt in order to get acceptable fits, the largest of 
which is 33 mA. Therefore, we conclude that the existing 
atomic absorption data are not sufficiently uncertain to mask the 
presence of this molecular signature. The search for molecular 
and solid features in X-ray spectra remains of great interest, 
both as tests of the atomic data and of models for the ISM. 

We believe that the present benchmark has produced an 
improved version of the warmabs model that can now be used for 
a more extensive and reliable study of oxygen K photoabsorption 
toward other X-ray sources. This will be the next step of the 
current study where the oxygen abundance variations in the ISM 
are certainly of interest. In a similar fashion, we also intend to 
evaluate the Ne and Fe edge regions for ongoing physical model 
refinement. 
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